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A  MECHANISM  FOR  IGNITION  OF  HIGH-TEMPERATURE  GASEOUS 
NITROMETHANE  -  THE  KEY  ROLE  OF  THE  NITRO  GROUP 
IN  CHEMICAL  EXPLOSIVES 

INTRODUCTION 

The  first  studies  of  nitromethane  pyrolysis  were  reported  by  Taylor 
and  Vesselovsky  [1]  and  Hirschlaff  and  Norrish  [2].  Considerably  later, 
plausible  initiation  steps  were  deduced  by  Cottrell,  Graham  and  Reid  [3]  and 
Hillenbrand  and  Kilpatrick  [4] .  The  first  accurate  measurement  of  the  rate  of 
unimolecular  decomposition  of  nitromethane  and  the  effects  of  pressure  was  done 
by  Glanzer  and  Troe  [S]  using  a  shock  tube.  This  work  established  that  C-N 
bond  rupture  is  the  only  plausible  initial  step.  Most  recently,  Perche, 

Tricot  and  Lucquin  [6]  made  the  first  attempt  to  construct  a  detailed 
chemical  kinetic  mechanism  for  the  decomposition  of  nitromethane.  Their 
scheme  was  applied  to  the  low-temperature  decomposition  of  nitromethane  vapor 
at  very  low  pressures. 

The  present  study  presents  a  chemical  reaction  mechanism  that  describes 
the  ignition  of  pure  gaseous  nitromethane  at  high  temperatures  and  identifies 
the  key  role  of  the  NO^  group  in  promoting  ignition.  This  work  was 
undertaken  as  a  first  step  in  understanding  the  chemistry  of  pure  solid 
explosives.  The  elementary  reaction  rate  constants  pertinent  to  the  decomposi¬ 
tion  of  nitromethane  have  either  been  collected  from  the  literature  or  estimated 
No  attempt  was  made  to  measure  any  new  rate  constant  or  to  confirm  a  previously 
measured  one.  Instead,  the  predictions  of  the  compiled  kinetic  mechanism 
were  tested  by  comparing  them  with  measured  induction  times.  Mass  and  infrared 
spectroscopy  were  used  to  identify  the  final  products.  Most  of  the  experiments 
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used  pure  nitrooethane  vapor,  although  some  used  mixtures  of  nitromethane 
diluted  In  argon  to  obtain  higher  temperatures  at  relatively  low  pressures. 
Temperatures  and  pressures  behind  the  reflected  shock  were  in  the  range  1000 
to  1600  K  and  1  to  10  atm. 

Traditionally,  shock  tube  experiments  with  100%  concentration  have  been 
avoided  because  of  the  large  departure  of  the  gasdynamic  configuration  from  one 
dimensionality.  In  such  cases,  calculating  the  rate  of  decomposition  requires 
complicated  corrections  to  the  effective  heating  time  and  temperature. 

However,  since  the  experiments  performed  here  measured  the  induction  time  only, 
the  pressure  at  the  center  of  the  end  plate  is  an  adequate  diagnostic 
provided  that  proper  care  is  taken  to  measure  the  incident  shock  velocity 
and  to  evaluate  the  temperature  behind  the  reflected  3hock.  The  details  of 
the  experimental  arrangement  and  diagnostics  are  given  below.  The  inherently 
one-dimensional  flow  character  at  the  center  of  the  end  plate  [7]  combined 
with  a  very  short  induction  time,  as  is  the  case  with  nitromethane  at  high 
temperatures,  allowed  the  use  of  a  one-dimensional  calculation  to  evaluate 
the  conditions  behind  the  reflected  shock. 

The  role  of  the  nitro  group,  NO^,  in  causing  rapid-exothermic  ignition 
was  tested  by  deleting  NO^  +  H  — >  OH  +  NO  from  the  set  of  reactions  describing 
the  decomposition  process.  Calculated  induction  times  and  final  temperatures 
were  then  compared  to  those  calculated  using  the  full  set  of  reaction  kinetics. 
In  all  cases,  the  induction  time  was  very  sensitive  to  the  above  reaction  rate. 
Below  we  describe  the  experiment,  and  discuss  the  chemical  kinetic  mechanism 
and  the  key  role  of  the  nitro  group  in  chemical  explosives. 


EXPERIMENT 


A  6.35  cm  internal  diameter  shock  tube  was  used  which  had  a  2.52  m 
driver  section  and  a  3.8  m  shock  section.  Helium  was  pumped  in  the  driver 
section  until  the  diaphragm  burst,  so  that  shocks  of  various  Mach  numbers  were 
achieved  by  changing  the  diaphragm  thickness.  Four  piezoelectric  gauges  located 
near  the  end  plate  detected  the  instant  of  shock  arrival  and  allowed  us  to 
evaluate  the  incident  shock  velocity  free  of  the  error  which  may  result  from 
its  attenuation.  A  CO  laser  probe  and  a  quartz  pressure  gauge  measured  the  NO 
absorption  and  the  pressure  evolution  at  a  fixed  location  4.5  cm  from  the  end 
wall.  A  second  pressure  gauge  measured  the  pressure  evolution  at  the  end  wall. 

After  all  of  the  wave  effects  in  the  shock  tube  died  out  and  mixing  was 
complete,  samples  of  the  products  were  extracted  for  later  analysis  by  mass 
spectrometry.  The  contents  of  the  shock  tube  were  then  pumped  through  a  liquid 
nitrogen  trap.  After  warming  up  the  contents  to  room  temperature,  an  infrared 
spectrometer  was  used  to  identify  the  condensible  products.  The  analysis  of 
collected  products  was  used  to  confirm  the  significant  species  formed  in  the 
decomposition  process.  No  attempt  was  made  to  measure  species  concentrations 
in  the  final  products.  Since  absorption  profiles  depend  on  temperature, 
the  measured  real-time  NO  absorption  gave  only  qualitative  evidence  of  the 
NO  behavior. 

Generally,  the  unshocked  gas  consisted  of  100%  spectro-grade 
nitromethane  at  pressures  up  to  a  few  torr.  In  a  few  cases,  argon  was  mixed 
with  nitromethane  so  that  the  mixture  could  be  shocked  to  high  temperatures 
without  increasing  the  driver  pressure  excessively  and  therefore  avoiding  damage 
to  the  piezoelectric  and  quartz  pressure  gauges.  However,  the  dilution  causes 
ignition  to  be  less  pronounced  on  the  pressure  trace  and  consequently  more 


difficult  to  detect 


The  conditions  behind  the  reflected  shock  were  calculated  from  the 
observed  incident  shock  velocity  assuming  frozen  composition  across  the  shock 
waves  and  equilibrium  of  the  internal  modes  of  the  nitromethane  molecules. 

The  effect  of  the  variation  of  the  heat  capacity  with  temperature  was 
included.  The  recorded  pressure  jumps  at  the  end  wall  were  found  to  compare 
favorably  with  those  calculated.  Temperatures  and  pressures  in  the  range  1000 
-  1600  K  and  1-10  atm  were  achieved.  The  reflected  shock  moves  into  a 
narrowing  channel  due  to  the  boundary  layer  formed  behind  the  incident  shock. 
Unlike  the  idealized  one-dimensional  problem,  in  this  case  the  gas  behind  the 
reflected  shock  does  not  come  to  zero  velocity.  Instead,  a  small  residual 
motion  slowly  compresses  the  gases  at  the  wall.  Because  of  the  very  short 
induction  times  of  nitromethane,  the  effect  of  this  compression  process  is 
ignored  and  a  constant  volume  adiabatic  process  starting  at  the  calculated 
temperature  and  pressure  behind  the  reflected  shock  is  assumed  for  the 
numerical  simulation.  However,  for  low-temperature  shocks,  when  the  heating 
due  to  the  subsequent  slow  compression  is  comparable  to  the  sudden  heating 
upon  crossing  the  shock,  the  induction  time  for  an  adiabatic  process  starting 
immediately  after  the  reflected  shock  can  be  significantly  longer  than  the 


measured  one 


REACTION  MECHANISM 


The  elementary  reactions  describing  the  decomposition  of  nitromethane 
at  high  temperatures  are  given  in  Table  I.  These  are  either  obtained  from  the 
literature  or  estimated.  Numerical  integration  [34]  of  the  rate  equations 
given  in  Table  I  gives  the  chemical  species  evolution  as  a  function  of  time. 
Tabulated  or  estimated  enthalpies  and  heat  capacities  [22a, 32]  of  the  chemical 
species  were  used  to  calculate  the  temperature  as  a  function  of  time.  Gas- 
dynamics  effects  were  not  calculated.  Instead,  an  adiabatic  constant  volume 
process  is  assumed. 

The  initiation  step  is  that  deduced  by  Glanzer  and  Troe  [S]  for  low 
pressures, 

CH3N02  +  M  -->  CH3  +  N02  +  M. 

The  major  overall  pathways  can  be  inferred  from  the  computations.  These  show 
that  the  reaction  proceeds  mainly  by  two  parallel  pathways,  as  shown  in 
Fig.  (1).  The  first  involves  H  abstraction  from  nitromethane  by  radicals 
to  form  CH2N02,  which  unlmolecularly  decomposes  into  formaldehyde.  The 
second  pathway  starts  with  the  methyl  plus  nitrogen  dioxide  reaction  to 
form  methoxy  radicals,  which  decompose  almost  instantaneously  into  formaldehyde, 
The  formaldehyde  from  both  pathways  is  then  attacked  by  radicals  to  produce 
formyl  radicals,  which  decompose  or  react  with  N02  to  produce  CO.  A  third 
minor  pathway  (not  shown  in  Fig.  (1))  involves  the  two  carbon  (  C2  )  species, 
whose  reactions  are  Included  in  Table  I. 

The  proposed  reaction  scheme  is  not  intended  to  be  all  inclusive.  Only 
those  reactions  (and  reverse  reactions)  having  rates  within  three  orders  of 
magnitude  of  the  fastest  rates  and  those  leading  to  observed  products  were 
retained  in  Table  I. 
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TABLE 


-  I  - 


Reaction  Rate  Constants,  k 


A  lVE/RT 


Reactants _ Products 


A _ a _ E _ Ref. 


Initiation 

CH3N02  +  M 

ch3  +  no2  +  m 

1.30(17) 

0.00 

42.000 

5 

Basic  Reactions 

ch3  +  N02 

CH30  +  NO 

1.30(13) 

0.00 

0.000 

5 

CH30  (M) 

CH.O  +  H  (M)a 

3.31(15) 

0.00 

27.500 

a 

ii 

c 

X  0.10615  ( 

[M]  T 

^0.65 

12 

N02  +  H 

NO  +  OH 

3.50(14) 

0.00 

1.500 

10 

CH30  +  OH 

ch2o  +  h2o 

3.20(13) 

0.00 

0.000 

8 

H2  +  OH 

H20  +  H 

2.20(13) 

0.00 

5.150 

15 

H2°  +  H 

H2  ♦  OH 

9.30(13) 

0.00 

20.400 

15 

ch3o  +  no2 

CH20  +  HONO 

4.00(11) 

0.00 

0.000 

9 

CH30  +  NO 

CH20  +  HNO 

3.2(12) 

0.00 

0.000 

9 

HONO  +  M 

NO  +  OH  +  M 

3.00(18) 

0.00 

46.700 

9 

NO  +  OH  +  M 

HONO  +  M 

7.91(15) 

0.00 

-2.200 

b 

HNO  +  M 

NO  +  H  +  M 

2.88(16) 

0.00 

48.800 

10 

NO  +  H  +  M 

HNO  +  M 

5.40(15) 

0.00 

-0.600 

10 

ch3no2  +  ch3 

CH2N02  +  CH4 

1.60(11) 

0.00 

10.800 

23 

CH3N02  +  H 

CH2N°2  ♦  H2 

6.30(13) 

0.00 

9.700 

30 

CH3N02  +  OH 

ch2no2  ♦  h2o 

1.85(12) 

0.00 

1.635 

17 

ch2no2  (M) 

CH20  +  NO  (M)8 

1.00(13) 

0.00 

36.000 

6 

Formvl  Droduction 

ch2o  +  ch3 

HCO  +  CH4 

3.10(10) 

0.00 

4.860 

31 

CH20  +  H 

HCO  +  H2 

1.26(13) 

0.00 

3.760 

25 

CH2C  +  OH 

HCO  +  HjO 

7.53(12) 

0.00 

0.175 

20 

(continued  on  next  page) 
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Products 


Products  formation 


HCO  +  M 

CO  +  H  +  M 

1.00(15) 

0.00 

14.700 

11 

CO  +  H  +  M 

HCO  +  M 

6.90(14) 

0.00 

-1.700 

16 

HCO  +  N02 

HONO  +  CO 

1.00(14) 

0.00 

0.000 

9 

HCO  +  NO 

HNO  +  CO 

2.00(11) 

0.50 

2.000 

8 

HCO  +  HNO 

CH20  +  NO 

3.20(13) 

0.00 

1.360 

8 

CO  +  OH 

C02  +  H  lo810k 

=  10.83  +  ■ 

3 . 94 ( —4 )  x 

T 

16 

CH .  +  OH 

4 

CH3  +  H2° 

1.45(12) 

0.00 

3.400 

19 

Two-Carbons  Grouo  Reactions 

ch3  +  ch3  (M) 

C2H6  (M)a 

3.16(13) 

0.00 

0.000 

27 

C2H6  +  CH3 

C2H5  +  CH4 

3.20(11) 

0.00 

10.800 

22 

C2H6  +  H 

C2H5  +  H2 

1.26(14) 

0.00 

9.700 

29 

C2H5  +  OH 

C2H5  +  H2° 

1.10(13) 

0.00 

2.440 

17 

C2H5  +  H  (M) 

CH3  +  CH3  (M)a 

3.71(13) 

0.00 

0.000 

26 

C2H5  +  M 

c2h4  +  H  +  m 

4.70(14) 

0.00 

26.600 

11 

C2H4  +  0H 

ch2o  +  ch3 

4.50(12) 

0.00 

0.200 

21 

C2H4  +  OH 

C2h3  +  H20 

7.33(12) 

0.00 

7.800 

17 

C2«3  +  M 

C2H2  +  H  +  M 

8.00(14) 

0.00 

31.500 

13 

Methanol  GrouD  Reactions 

CH3  +  OH  (M) 

CHjOH  (M)a 

8.00(12) 

0.00 

0.000 

11 

ch3oh  +  ch3 

ch2oh  +  ch4 

3.15(10) 

0.00 

6.360 

31 

ch3oh  +  H 

ch2oh  +  H2 

1.50(13) 

0.00 

5.260 

24 

CH30H  +  OH 

CH20H  +  H20 

4.23(12) 

0.00 

0.850 

18 

ch2oh  (M) 

CH30  (M)a 

1.00(13) 

0.00 

39.500 

28 

CH20H  +  H 

CH20  +h2 

7.50(13) 

0.00 

0.000 

24 

(continued 

on  next 

page) 
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TABLE  -  I  -  (continued) 


Reactants 

Products 

A 

n 

E 

Ref, 

Nitrosomethane  GrouD  Reactions 

CH3  +  NO  (M> 

CH3NO  (M)a 

4.00(12) 

0.00 

0.000 

6 

ch3no  (M) 

CH3  +  NO  (M)a 

7.00(13) 

0.00 

38.000 

6 

ch3no  +  ch3 

ch2no  +  ch4 

1.60(11) 

0.00 

10.800 

23 

CH3NO  +  H 

ch2no  +  h2 

6.30(13) 

0.00 

9.700 

30 

CHjNO  +  CH 

CH2NO  +  H20 

1.85(12) 

0.00 

1.635 

17 

ch2no  (M) 

HCN  +  OH  (M)a 

3.98(11) 

0.00 

33.000 

14 

(a)  Rate  constants  of  reactions  with  ' (M) '  have  a  pressure  fall-off  curve. 
Units  of  rate  constant  are  indicated  by  number  of  reactant  species,  where 

the  parenthetical  (M)  does  not  count.  For  example,  the  rate  constant 

3  -1 

for  CH  NO  +  M  has  units  cm  /mol  s,  that  for  CH.O  (M)  has  units  s  , 

*  2  6  2* 
and  that  for  NO  +  OH  +  M  has  units  cm  /mol  s.  Activation  energy  is 

given  in  Kcal/mol. 

(b)  Reverse  reaction  rate  evaluated  from  preceeding  forward  reaction 
rate  and  equilibrium  thermochemistry. 
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Fig.  3  —  Most  significant  reaction  rates  in  the  decomposition  of  100%  nitromethane,  calculated  for  an 
adiabatic  constant  volume  process  (a)  starting  at  1108  K  and  3.92  atm,  and  (b)  starting  at  1494  K  and 
5.73  atm. 
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(c)  Pressure  trace  at  end  wall;  7.7  torr  100%  nitromethane;  incident  shock  velocity:  0.122  cm/^ts; 
calculated  temperature  and  pressure  behind  reflected  shock:  1227  K  and  4.82  atm;  measured  induction 
time:  10-12  (is. 
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(d)  Pressure  trace  at  end  wall;  10.31  torr  100%  nitromethane;  incident  shock  velocity:  0.125  cm//is; 
calculated  temperature  and  pressure  behind  reflected  shock:  1263  K  and  6.96  atm;  measured  induction 
time:  5  n$. 

Fig.  2  vCont’d)  —  Pressure  and  NO  absorption  traces  from  the  shock  tube  experiment.  Measurements 
taken  either  at,  or  4.5  cm  from,  the  end  wall.  Numbers  at  the  top  refer  to  oscilloscope  gain  per  division 
for  pressure  and  NO  traces,  respectively.  Third  number  is  the  time  increment  per  division.  Observed 
oscillations  in  pressure  occuring  after  the  first  spike,  denoting  ignition,  should  not  affect  induction  time 
increments. 


(a)  Pressure  (lower)  and  NO  (upper)  traces  4.5  cm  from  end  wall;  mixture:  5.0  torr  nitromethane, 
14.5  torr  argon;  incident  shock  velocity:  0.086  cm/(is;  calculated  temperature  and  pressure  behind 
reflected  shock:  1145  K  and  1.69  atm;  measured  induction  time:  70-80  ft s. 


WCm* 


(b)  Pressure  (lower)  trace  at  end  wall;  NO  (upper)  trace  4.5  cm  from  end  wall;  18.19  torr  100% 
nitromethane;  incident  shock  velocity:  0.112  cm  /(is\  calculated  temperature  and  pressure  behind 
reflected  shock:  1108  K  and  8.8  atm;  measured  induction  time:  25-30  (is. 

Fig.  2  —  Pressure  and  NO  absorption  traces  from  the  shock  tube  experiments.  Measurements  taken 
either  at,  or  4.5  cm  from,  the  end  wall.  Numbers  at  the  top  refer  to  oscilloscope  gain  per  division  for 
pressure  and  NO  traces,  respectively.  Third  number  is  the  time  increment  per  division.  Observed 
oscillations  in  pressure  occuring  after  the  first  spike,  denoting  ignition,  should  not  affect  induction  time 
measurements. 


pathways  for  the  decomposition  of  gaseous  nitromethane 


CONCLUSION 


The  decomposition  of  CHjN02  proceeds  by  two  major  parallel  pathways, 
through  methoxy  radicals  and  through  CH2N02  radicals.  Both  pathways  produce 
CH20  which,  when  attacked  by  OH  and  H  radicals,  yield  HCO  radicals  that  carry 
the  reaction  towards  completion.  As  long  as  unreacted  nitromethane  remains, 
its  unimolecular  (C-N  bond  rupture)  reaction  serves  as  an  energy  sink, 
preventing  ignition.  As  a  result,  although  the  radical  concentrations  reach 
their  maximum  level  early  in  the  overall  process,  ignition  does  not  occur 
until  all  nitromethane  is  consumed.  The  energy  released  then  by  the  radical 
reactions  causes  rapid  ignition. 

The  computations  show  that  the  reaction 

N02  +  H  — >  OH  +  NO 

plays  a  key  role  in  accelerating  ignition  in  nitromethane  because  it  controls 
the  OH  concentration,  due  to  the  absence  of  other  efficient  routes  to  OH. 

This  effect  is  expected  in  other  explosives  containing  C,  H  and  the  nitro 
or  nitrate  [9]  groups,  to  the  extent  that  other  efficient  OH  formation  routes 
are  lacking. 

Ve  have  compared  the  predictions  of  the  chemical  kinetic  scheme  to  the 
shock  tube  experiments.  The  mechanism  proposed  gives  the  correct  induction 
times  in  the  pressure  and  temperature  ranges  considered. 
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Calculated  Effect  of  k(H  +  NO  — >  OH  +  NO)  on  Induction  Time 
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ignition  did  not  occur 


The  H  atom  concentration  was  approximately  doubled  by  setting  the 
rate  of  H  +  N02  — >  OH  +  NO  equal  to  zero.  The  resulting  increase  in  the  rates 
of  the  remaining  H  atom  reactions  did  not  shorten  the  calculated  ignition  time. 
This  fact,  together  with  the  rates  plotted  in  Figs.  (3a)  and  (3b),  shows  that 
OH  is  the  key  species  consuming  nitromethane  in  the  radical  chain  mechanism. 

It  is  also  interesting  to  note  that  the  heat  released  by  methyl  radical  re¬ 
combination  did  not  promote  ignition,  unless  the  radical  chain  mechanism 
initiated  by  H  +  NO^  — >  OH  +  NO  was  included. 

The  effect  of  varying  other  selected  rate  constants  was  also  tested, 
however  systematic  sensitivity  analyses  were  not  done.  In  general,  the  reactions 
with  the  largest  rates  will  control  the  decomposition  process,  and  these  are 
shown  in  Figs.  (3a)  and  (3b).  These  reactions  have  well  established,  measured 
rate  constants,  with  a  few  exceptions.  The  estimated  Arrhenius  parameters  of 
Perche  et.  al.  for  unimolecular  decomposition  of  CH^NO^  [6]  could  give  a  rate 
that  is  in  error  by  a  factor  as  large  as  ten.  Nevertheless,  the  lifetime  of  this 
radical  is  short  with  respect  to  unimolecular  decay,  and  this  decay  path  is  the 
dominant  one  for  its  consumption.  We  have  made  estimates,  which  should  be  cor¬ 
rect  to  within  a  factor  of  three,  for  the  rate  constants  of  the  hydrogen 
abstraction  reactions  of  H  and  OH  with  nitromethane  [17,30].  Confidence 
in  these  estimates  is  based  upon  agreement  between  predicted  and  measured  rate 
constants  for  H  and  OH  abstracting  hydrogen  from  C-H  bonds  in  other 
molecules  [17,19,22].  Finally,  methyl  radical  recombination  to  form  ethane, 
and  methoxy  radical  unimolecular  decomposition,  have  pressure  dependent  rate 
constants.  Appropriate  pressure  adjusted  values  have  been  used  [12,33]. 


ROLE  OF  NITRO  GROUP 


The  sensitivity  of  the  calculated  induction  time  to  the  reaction. 


N02  +  H  ~>  OH  +  NO 

was  tested  by  deleting  it  from  the  set  of  reactions  describing  nitromethane 
decomposition.  The  results,  for  a  representative  subset  of  the  initial  conditions. 


are  given  in  Table  III.  Deleting  this  reaction  drastically  increased  the 
induction  time  at  lower  temperatures,  the  difference  becoming  smaller  as  the 
temperature  increased.  For  example,  while  nitromethane  ignites  in  35  -  40  (is 


at  1108  K  and  3.92  atm,  no  sign  of  ignition  was  observed  after  200  (is  when 
the  above  reaction  was  deleted.  At  1263  K  and  6.96  atm,  ignition  had  not 
occured  after  20  ns  compared  to  2  -3  (is  for  the  full  set  of  reactions.  The^ 
Induction  time  only  triples  from  0.4  to  1.2  (is  at  1494  K  and  5.73  atm.  However, 
although  ignition  occured,  the  overall  reaction  system  remained  globally 
endothermic.  The  final  calculated  temperature,  1430  K,  is  less  than  the  starting 


temperature.  1494  K,  and  much  smaller  than  the  final  temperature  calculated  with 
the  full  set  of  reactions,  1840  K. 


The  rapid  ignition  of  nitromethane  is  attributed  to  the  abundance  of  N02 
formed  early  in  the  decomposition  process,  specifically  to  the  rapid  reaction 
of  N02  with  H  radicals  producing  large  quantities  of  OH.  The  subsequent 
chain-carrying,  H-abstraction  reactions  of  OH  with  nitromethane  and  formalde¬ 
hyde  are  fast  and  exothermic  and  contribute  to  the  release  of  energy.  The 
reaction  of  OH  with  nitromethane  accelerates  the  consumption  of  starting 
material,  nitromethane.  The  formaldehyde  reaction  with  OH  leads  to  the 
regeneration  of  H  which  restarts  the  chain. 
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TABLE  -  II  - 

Calculated  and  Observed  Induction  Times 


Data 

Pressure 

Temperature 

Induction 

Time  (ps) 

Point  *  CH3N02 

(atm) 

(K) 

Calculated 

Measured 

100%  Nitromethane 

1  100% 

4.00 

1017 

>240a 

no  ignition 

2 

4.00 

1035 

>150a 

no  ignition 

3 

5.78 

1030 

>150a 

80 

4 

3.92 

1108 

35-40 

35-40 

5 

4.00 

1120 

30-35 

35 

6 

5.20 

1142 

15-20 

20  and  30-: 

7 

4.96 

1120 

25-30 

30 

8 

3.76 

1155 

15-20 

25-30 

9 

8.80 

1108 

20-25 

25-30 

10 

7.25 

1131 

15-20 

12  and  25° 

11 

8.17 

1180 

5-10 

7 

12 

4.82 

1227 

5-6 

10-12 

13 

6.85 

1202 

6 

10-12 

14 

6.96 

1263 

2-3 

5 

15 

6.56 

1283 

2.2 

4 

16 

6.17 

1376 

0.8-1 

2_3d 

17  5.73 

Mixtures  of  Nitromethane  and  Arson 

1494 

0.4 

<  ld 

1  26 

1.69 

1145 

60-70 

70-80 

2  64 

3.94 

1341 

2-3 

3-5 

3  49 

2.18 

1440 

2-3 

4 

4  42 

4.52 

1542 

0.5-1 

1-2 

a  Temperature  Increased  slowly.  No  sudden  rise  In  temperature  was 
observed  up  to  listed  time, 
b  No  ignition  spike  was  observed  on  pressure  trace, 
c  Two  Ignition  spikes  were  observed  on  pressure  trace.  The  second 
was  higher  in  amplitude  and  sharper, 
d  Ignition  spike  undetectable  from  reflected  shock  jump  on  the 
oscilloscope  pressure  trace. 
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Although  most  of  the  radicals  reach  their  maximum  concentration  early 
in  the  reaction  process,  ignition  does  not  occur  until  all  of  the 
nitrome thane  is  consumed.  As  long  as  nitrome thane  is  still  present,  the 
overall  reaction  system  stays  endothermic  because  the  radical  reactions  do 
not  release  enough  energy  to  compensate  for  the  energy  consumed  in  breaking 
the  C-N  bond.  When  the  nitromethane  is  depleted,  the  system  suddenly 
becomes  exothermic  and  ignition  occurs.  In  this  mechanism,  all  reaction 
rates  drop  sharply  at  ignition  without  the  preceding  sudden  increase 
observed  in  most  systems.  This  behavior  is  due  to  the  sudden  decrease  in  the 
concentrations  of  reactive  species  which  have  rapidly  reacted  to  form  stable 
products,  without  any  nitromethane  present  to  replenish  them. 

Figure  (4a)  shows  the  calculated  time  evolution  of  the  intermediate 
species  for  the  same  initial  conditions  as  in  Fig.  (2d).  The  methoxy  and 
hydrogen  radicals,  which  are  consumed  almost  as  fast  as  they  are  formed,  have 
a  very  low,  nearly  steady-state  concentration.  We  note  an  abundance  of  N02, 
which,  as  explained  below,  accelerates  the  ignition  process.  Formaldehyde  has 
the  highest  concentration  of  any  intermediate,  and  its  depletion  begins  late 
in  the  induction  period.  Figure  (4b)  shows  the  concentrations  of  stable 
products  steadily  increasing  until  the  end  of  the  induction  time,  2.5  ps, 
when  they  level  off.  In  our  calculations,  this  time  coincides  consistently 
with  that  of  consumption  of  all  the  nitromethane.  We  note  the  same  behavior 
of  the  NO  absorption  trace  in  Fig.  (2a). 


shock  and  the  Ignition  spike  on  the  oscilloscope  pressure  trace.  By  fitting 
the  calculated  induction  times  for  100%  nitrometbane  to  an  Arrhenius  form  we 
obtain 

x^  =  Aiexp(Ei/RT) 

where 

Ai  -  6.13  x  10"7  |is 
E±  *  38.6  Kcal/mol. 

For  temperatures  above  1600  K  the  calculations  should  follow  this  same 
expression.  We  believe  that  the  proposed  mechanism  should  correctly  predict 
ignition  behavior  up  to  2000  K. 

The  calculated  rates  of  the  most  significant  reactions  for  the  initial 
conditions,  1108  K,  3.92  atm  and  1494  K,  S.73  atm,  corresponding  to  points 
4  and  17  in  Table  II,  are  shown  in  Figs.  (3a)  and  (3b),  respectively.  The 
same  sequence  of  events  is  observed  in  the  computations  at  both  temperatures: 

*  C-N  bond  in  CH3N02  breaks  giving  CHj  +  N02 

*  CH3  radicals  recombine  producing  CjH^ 

*  CH3  recombines  with  N02  producing  CHjO,  which  immediately  decomposes  giving 
H  radicals  +  CHjO,  and  NO 

*  H  radicals  attack  CH3N02  producing  CH2N02  radicals 

*  H  radicals  react  with  N02  producing  OH  radicals  and  NO 

*  OH  radicals  attack  CH3N02  producing  CHjNOj  radicals 

*  CH2N02  radicals  unimolecularly  decompose  giving  more  CH20 

*  OH  (and  H)  radicals  attack  CH20  giving  HC0 

*  HC0  radicals  decompose  or  react  with  N02  to  give  CO. 


RESULTS 


Figure  (2a)  illustrates  typical  pressure  (bottom)  and  NO-absorption 
(top)  oscilloscope  traces  at  the  fixed  location  4.5  cm  from  the  end  wall,  for 
a  mixture  of  26%  nitromethane  in  argon.  The  first  two  numbers  at  the  top  of  the 
figure  refer  to  the  oscilloscope  gain  (pressure,  NO-absorption)  while  the  third 
one  refers  to  the  time  increment  per  division.  The  sudden  jump  in  the  NO 
absorption  trace  marks  the  arrival  time  of  the  reflected  shock.  At  70-80  (is, 
after  the  pressure  jump  caused  by  the  reflected  shock,  there  is  another 
pressure  jump  that  signals  the  end  of  the  induction  period.  This  coincides  with 
the  NO  trace  suddenly  leveling  off.  Figures  (2b-2d)  illustrate  typical  pressure 
traces  recorded  at  the  end  wall  for  pure  nitromethane.  The  pressure  jump  in 
each  record  designates  the  reflected  shock  while  the  pressure  spike  indicates 
ignition.  The  latter  is  more  pronounced  for  higher  reflected  shock  pressures. 

The  induction  times  are  25-30  (is,  10-12  ns,  and  5  |is,  respectively.  The  slow 
pressure  rise  following  the  reflected  shock  in  Fig.  2c  is  due  to  the  residual 
particle  motion  towards  the  end  wall,  typical  when  the  reflected  shock  sweeps 
through  a  channel  narrowed  by  the  growth  of  the  boundary  layer  on  the  side  walls. 
The  final  products,  HjO,  CO,  CO^,  NO,  N02,  HCN,  C^H^,  C^H^,  were  consistently 
observed . 

Measured  induction  times  and  those  calculated  by  Integrating 
the  rate  equations  for  an  adiabatic  constant  volume  process  are  listed  in 
Table  II  for  different  initial  conditions  ranging  from  1000  to  1600  K  and  1  to 
10  atm.  Most  of  the  data  listed  are  for  100%  nitromethane.  The  results  for 
temperatures  above  1600  K  (behind  the  reflected  shock)  were  disregarded  because 
the  incident  shock  is  then  strong  enough  to  cause  a  significant  decomposition  of 
the  starting  material.  Moreover,  for  temperatures  above  1600  K  the  induction 
time  is  so  short  that  it  was  impossible  to  distinguish  between  the  reflected 


the  estimation  of  Arrhenius  parameters  for  the  isomerization  of  alkoxy 


radicals  to  hydroxyalkyl  radicals,  a  structurally  analogous  change.  We 
have  assumed  a  strain  energy  of  26  kcal/mole  for  the  4-membered  ring 
transition  state  in  this  case. 

15.  Baulch.  D.,  Drysdale,  D.,  Horne,  D.,  and  Lloyd,  A.,  Evaluated  Kinetic 
Data  for  High  Temperature  Reactions,  Vol.  I,  The  University  of  Leeds,  1972. 

16.  Baulch,  D.,  Drysdale,  D.,  Duxbury,  J.,  and  Grant,  S.,  Evaluated  Kinetic 
Data  for  High  Temperature  Reactions,  Vol.  Ill,  The  University  of  Leeds,  1976 

17.  (a)  Greiner,  N.  R.,  J.  Chem.  Phys.,  53,  1070  (1970).  This  reference 
reports  several  measured  H-abstraction  rate  constants  for  OH  and  derived 
expressions  for  OH  reaction  rates  with  primary,  secondary,  and  tertiary 
R-H  bonds.  The  following  references,  17b  -  17e,  describe  our  use  of  this 
data. 

(b)  For  CH^NO  and  CHgN02,  a  primary  C-H  bond  energy  was  assumed,  and  A 
and  E  were  calculated  from  Greiner's  expression. 

(c)  Greiner's  measurement  of  the  rate  constant  for  ethane  was  used. 
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